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Abstract. A global database of approximately 1.7 million observations 
of the partial pressure of carbon dioxide in surface ocean waters ( pCOs ) col- 
lected between 1970 and 2003 is used to estimate its spatial autocorrelation 
structure. The patterns of the lag distance where the autocorrelation exceeds 
0.8 is similar to patterns in the spatial distribution of the first baroclinic Rossby 
radius of deformation indicating that ocean circulation processes play a sig- 
nificant role in determining the spatial variability of pCO s . For example, the 
global maximum of the distance at which autocorrelations exceed 0.8 aver- 
ages about 140 km in the equatorial Pacific. Also, the lag distance at which 
the autocorrelation exceed 0.8 is greater in the vicinity of the Gulf Stream 
than it is near the Kuroshio, approximately 50km near the Gulf Stream as 
opposed to 20 km near the Kuroshio. Separate calculations for times when 
the sun is north and south of the equator revealed no obvious seasonal de- 
pendence of the spatial autocorrelation scales. The pCO s measurements at 
Ocean Weather Station (OWS) “P” in the eastern subarctic Pacific (50°N, 
145°W) is the only fixed location where an uninterrupted time series of suf- 
ficient length exists to calculate a meaningful temporal autocorrelation func- 
tion for lags greater than a few days. The estimated temporal autocorrela- 
tion function at OWS “P” is highly variable. A spectral analysis of the longest 
four pCOs time series indicates a high level of variability occurring over pe- 
riods from the atmospheric synoptic to the maximum length of the time se- 
ries, in this case 42 days. It is likely that a relative peak in variability with 
a period of 3-6 days is related to atmospheric synoptic period variability and 
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ocean mixing events due to wind stirring. However, the short length of avail- 
able time series makes identifying temporal relationships between pCC>2 and 
atmospheric or ocean processes problematic. 
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1. INTRODUCTION 

The partial pressure of C 0 2 (pCO s ) in seawater is a vapor pressure of C 0 2 , and hence 
governs the magnitude of the C0 2 transfer flux across the sea-air interface when multiplied 
with sea-air gas transfer coefficient. It is a sensitive function of temperature doubling 
with every 16 °C (Takahashi et al., 1993). It is also a sensitive function of the total 
concentration of C0 2 (TC0 2 ) species dissolved in seawater 

TC0 2 = [C0 2 ] aqueous + H C O 3 + CO 3 

that depends on the net biological community production, the rate of upwelling of subsur- 
face waters rich in C0 2 , and the air-sea C0 2 flux. The local sensitivity may be expressed 
in terms of the Revelle factor, q\™tco 2 , which varies from 8 with lower TC0 2 concentra- 
tions in tropical waters to 15 with higher concentrations in polar waters. In the surface 
mixed layer, the effect of warming on pCOg is counteracted by lower TC0 2 caused by 
photosynthetic fixation of C0 2 , as often seen during spring bloom periods; and the effect 
of cooling is counteracted by increasing TC0 2 usually caused by upwelling of subsurface 
waters rich in respired C0 2 . Consider the following example. A parcel of polar ocean 
water at -1.9 °C warmed to an equatorial temperature of 30°C without changes in TC0 2 
and other chemicals, increases its pCOs by a factor of 4. If the typical nitrate in this 
water (~35 /xmol/kg) is completely utilized by biological growth with the Redfield N/C 
ratio of 16/106, then the TC0 2 in the same water would decrease from 2150 pmol/kg to 
1920 pmol/kg. Thus, there is a decrease in pCO s by a factor of (1920/2150) 10 r* 33 % 
(using a Revelle factor of 10). This example illustrates that, over the global oceans, the 
effect of change in temperature is roughly compensated by changes in TC0 2 , and that the 
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time-space variation in surface seawater pC0 2 is dictated primarily by interactions be- 
tween the effects of temperature, net biological production and the deep water upwelling 
rate. Using climatological mean data, Takahashi et al. (2002) have observed that, in 
subtropical gyres, seasonal changes in pC0 2 are due primarily to temperature changes, 
whereas those in sub-polar and polar waters are due primarily to TC0 2 changes caused 
by winter upwelling of subsurface waters and spring plankton blooms in spring. However, 
finer scale time-space variability of surface water pC0 2 has not been investigated using 
semi-continuous underway surface water pC0 2 measurements which have been made in 
increasing numbers during recent years. 

In this study, the global distribution of the isotropic spatial autocorrelation structure 
of available pC0 2 underway measurements is calculated. The motivation is to identify 
the scales of variability that are resolved in the pC0 2 data, and attempt to relate some 
of that variability to known physical processes. In addition, time series data available 
from OWS “P” is similarly analyzed to provide an example of temporal autocorrelation 
structure and variability. It is hoped that a quantified estimate of autocorrelation scales 
will be useful to investigators, particularly future data assimilation efforts, for decisions 
on how best to use the 1.7 million measurements in choosing averaging and sub-sampling 
schemes for their work. 

2. DATA 

For the global studies, most of the pC0 2 measurements in surface mixed layer wa- 
ters (hereafter referred to simply as pC0 2 ) are mostly comprised of underway shipboard 
measurements and a much smaller number of measurements for discrete water samples 
collected at hydrographic stations. The underway data used for both Takahashi et al. 
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(1997, 2002) studies was temporally averaged to 12-hourly intervals; the averaging reduc- 
ing the number of total observations as underway measurements in the complete data 
set are sometimes taken as frequently as every 2 minutes, but usually several times per 
hour. The 12-hourly (approximately 100 km) averaging period was chosen ostensibly to 
maintain the larger scale gradient structure of pCOg in the total data set. 

While the raw number of pCOz measurements (currently approaching 1.7 million in 
some data sets) may seem adequate, or perhaps even large, for global studies, the number 
of actual data values used in past studies is much lower because of the averaging of 
the underway data as noted above. Even the choice of averaging period/distance is, 
in some ways, an ad hoc choice depending on the scales of interest and the accuracy 
of underway data versus station point data. Indeed, the relative paucity of accurate 
pCOg measurements has lead some investigators to try to infer pCOg values from other 
properties of the seawater, such as temperature. For example, Boutin et al. (1999) and 
Feely et al. (2004) found strong co- variability between pCO s and sea surface temperature 
in the equatorial Pacific, and Loukos et al. (2000) were able to show a useful predictive 
skill of temperature and salinity for TCO 2 in the equatorial Pacific. 

The pC0 2 data used here is an expanded version of the data used in Takahashi et al. - 
(2002), including about 760,000 new observations bringing the total to approximately 1.7 
million. The data was collected between 1970 and 2003, and its geographic and rough 
temporal distributions are shown in Fig. 1. Note, for the most part, the entire Indian 
and South Pacific Oceans are only sparsely sampled prior to 1990. Much of the sampling 
in the South Atlantic occurred between 1985 and 1989. Prior to 1985, the only sub basin 
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that is sampled more densely is the North Atlantic. A caution on the results presented 
here is that the results apply to sampling that primarily took place post 1984. 

All autocorrelation calculations are performed on data interpolated to regularly spaced 
segments with a distance of 2 km between successive data points. Interpolating pC0 2 data 
across gaps larger than 20 km is not done. For gaps larger than 20 km, a new segment 
is defined. Also, interpolation across a gap in time of longer than 10 days is not allowed, 
and a new segment would be defined. This latter condition was very rarely met. Average 
autocorrelation functions are calculated globally on a 10° x 10° grid. The location of the 
segment is taken to be the average latitude and longitude of the measurement locations. 
Segments shorter than 60 km were not included in the calculations, and no single segment 
was allowed to exceed 1000 km, approximately the length of the grid box. For segments 
longer than 1000 km, a group of subsegments, broken every 1000 km, was used for the 
calculations. Thus, the autocorrelations describe variability on approximate spatial scales 
less than the grid box resolution. 

The global distribution of segments is not homogeneous, and the number of segments 
in each grid box is presented in Fig. 2. Areas of densest coverage include the tropical 
Pacific, the North Pacific near the coast of Japan and north of 45°N, the northwest Indian 
Ocean, the Caribbean and specific locales in the Southern Ocean. In many places, the 
correlations presented are the result from only 1 or 2 segments, such as is the case for much 
of the South Atlantic and Pacific Oceans. In those areas, little confidence can be placed 
in the autocorrelation structure, but they are nevertheless presented for completeness. 

In addition to the global 10° grid, seven ocean-relevant areal averages of autocorrelation 
functions are presented. Those areas include: the Kuroshio subtropical gyre, 135°E-180°E, 
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25°N-40°N; eastern equatorial Pacific, 180°E-80°W, 10°S-10°N; northeast Pacific, 150°W- 
120°W, 30°N-50°N; East Australia gyre, 150°E-180°E, 40°S-10°S; Gulf Stream subtropical 
gyre, 80°W-50°W, 25°N-40°N; north Atlantic subpolar gyre, 60°W-10°W, 40°N-60°N; and 
the tropical Indian, 30°E-100°E, 20°S-10°N. 

3. SPATIAL AUTOCORRELATION STRUCTURE 

The lag autocorrelation function, AC(lag), is calculated as follows: 

a/~i( 1 \ _ Si=0 9 ( x i ~ x ){ x i+lag ~ x ) 

AC{la9) ~ 

where a: is the measurement value and an overbar represents the average value. The spa- 
tial distance of the lag is simply the product of the index lag times the grid resolution 
between consecutive data points, in this case 2 km. For this study, we present results for 
an isotropic autocorrelation only, understanding that there are locations where the auto- 
correlation function is almost certainly not isotropic, near 40° - 45°N in the North Pacific, 
e.g. Murphy et al. (2001) performed identical calculations for 24 basin long repeat tracks 
of pCOg measurements made over a 48 months period along nearly great circle routes 
in the North Pacific. That study concentrated on the large scale spatial decorrelation 
estimates, and found that pCOg data became completely decorrelated (AC = 0) at about 
1.5° longitude (105 km along the E-W direction) in the Bering Sea and 3.5° longitude (250 
km along the NW- SE direction) in the Gulf of Alaska. In this study, it is the smaller 
scale variability that is of interest, and many data segments are not of sufficient length 
that would allow AC to approach 0. 

Rather than presenting an estimate of the spatial autocorrelation function for each grid 
box, the global structure is presented as the lag distance at which the autocorrelation 
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estimate drops below a specific value. The first map (Fig. 3a) indicates the distance over 

which the autocorrelation is termed high, here, greater or equal to 0.8. The second map 
is the distance over which the autocorrelation is termed elevated, 0.6. These definitions of 
high and elevated are arbitrary, but nevertheless convenient for discussion purposes. The 
statistical significance of each autocorrelation is a function proportional to the number of 
data points that enter into the calculation of AC and varies geographically in this study 
(see Fig. 2). The globally recognizable feature is that the lag distance is shorter in high 
latitude areas and increases toward the equator, consistent with the dependence of the 
Rossby radius of deformation on the inverse of the local Coriolis acceleration. 

The largest area where high autocorrelations extend to the greatest distance is in the 
immediate vicinity of the equatorial Pacific where that distance is routinely greater than 
80 km, ranging from 100 km to 250 km in thermocline waters. There is some indication 
this “high correlation distance” is greater north of the equator than it is to the south in the 
Pacific Ocean. In the Atlantic, the high correlation distance is greater north of the equator 
than it is directly on the equator. In the Indian Ocean, the high correlation distance is 
greater south of the equator than it is to the north of the equator. The asymmetry in the 
Indian Ocean may be due to the long annual Rossby waves that dominate much of the 
variability in the southern Indian Ocean but does not in the north. 

In the northern hemisphere subtropical gyres, the distance is greater (about 50 km) in 
the Atlantic compared to the Pacific where the distance is below 30 km. This distribution 
would be consistent with the eddy transfer scales in these gyres as the internal deformation 
radius of the subtropical gyre in the Atlantic is about 45 km, and that deformation radius 
is about 25 km in the Pacific (see Emery et al, 1983 for estimates of deformation radii 
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based on local vertical Brunt- Vaissala profiles). At high latitudes, the high correlation 
distance is 20 km and shorter, consistent again with the smaller deformation radii at these 
latitudes. Thus, in general, the distance for autocorrelations greater than 0.8 is consistent 
with scales set by the physical properties of the ocean flow in many areas. 

Not surprisingly, there are aspects of the elevated autocorrelation distance (Fig. 3b) 
that are similar to the patterns of high correlation distance. For example, the longest 
distance for elevated autocorrelations is in the equatorial Pacific. Also, there are longer 
distances for elevated autocorrelation in the Gulf Stream gyre compared to the Kuroshio 
subtropical gyre are as in the high correlation distance pattern. However, there are areas 
in the high latitudes that have elevated autocorrelation distance that are comparable to 
the distances in the subtropics. This is somewhat evident in both the North Atlantic and 
North Pacific, indicating the autocorrelation function does not have a steep drop between 
0.8 and 0.6 in those areas. Also note that the distances with elevated correlation are 
longer in the Gulf of Alaska compared to those east of the dateline at the same latitude 
consistent with the calculations of Murphy et al. (2001). 

It is worth noting that even though the Murphy et al. (2001) study and this study used 
differing resolutions for the input data (10 km vs 2 km for the present study), both studies 
seemingly arrive at the same lag distance over which the autocorrelation drops below 0.8 
in the northern North Pacific, about 20-40 km. However, the same zonal asymmetry 
appearing in both studies is a little surprising. The data used by Murphy et al. (2001) for 
the northernmost tracks between 170°E and 165°W were restricted within the Bering Sea, 
and the pC0 2 variability may likely reflect local turbulence generated by the Aleutian 
Island chain. In contrast, the variability of the northern North Pacific in our study is 
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represented mainly by data south of the island chain in the open North Pacific. Here, the 
shortened autocorrelation scales seem more likely tied to the variability of the Kuroshio 
Extension and Oyashio as opposed to local island effects that may have been important 
in the analysis of Murphy et al. (2001). 

Around the Antarctic Circumpolar Current, the situation is more complex and changes 

1 

in the elevated autocorrelation distance often varies greatly in adjacent grid boxes. The 
variability in this region is controlled by a strongly meandering zonal flow. The spatial 
scale of the local variability is much smaller than the 10° x 10° box used in this study. 
Thus the Antarctic Circumpolar Current and its variability are very much under-sampled, 
and this may explain the strong box to box variability in the autocorrelation calculations 
presented here. 

The same calculations were repeated, except that average autocorrelation structures 
were calculated for two seasons, March-August and September-February, when the sun is 
at maximum zenith in the northern and southern hemisphere, respectively. Except for a 
slight tendency for there to be longer autocorrelation distances in the subtropical gyres 
when the sun is at zenith in the local hemisphere, there is very little seasonal change with 
the calculated autocorrelation structure. The change in the subtropical gyres was not sig- 
nificant enough to warrant showing the figures for the seasonal changes. The rapid change 
in autocorrelation structure over relatively short distances near the Antarctic Circumpo- 
lar Current did not change with season. Thus, biology seems an unlikely candidate for a 
single controlling factor of this spatial behavior. This is consistent with the small seasonal 
amplitudes of pGO s observed in the 40°S-60°S zone suggesting that pCOs changes due 
to SST are largely compensated by changes in TCO 2 due to the biological utilization and 
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upwelling (Takahashi et al., 2002). The highly variable current and wind forcing make it 
difficult to assign any one reason for such variability in this area of strong horizontal and 
vertical mixing. 

Fig. 4 is the estimated autocorrelation structure for seven specific geographic areas 
defined in the DATA section of this manuscript. Along with the estimated structure, a 
least squares best fit exponential curve to the average autocorrelation function is given to 
provide an estimate of an exponential decay scale of this structure. As indicated by the 
10° x 10° estimates, the equatorial Pacific has the longest exponential decay scale of about 
139 km. Also, the Gulf Stream decorrelation scale is nearly double than that occurs near 
the Kuroshio. An unexpected result is the relatively long 92 km e-folding scale in the 
North Atlantic subpolar gyre. Note, however, that the drop off for the highest correlation 
scales is not represented well by the best fit line in this area and a simple exponential decay 
would be an obviously poor choice for a functional representation of the autocorrelation 
function. This is also the case in the East Australia current region. Nevertheless, the 
exponential estimate provides consistent, if overly simplistic, comparison for the regions. 

4. TEMPORAL AUTOCORRELATION STRUCTURE 

Within this pCOs database, there are but a few time series of sufficient length to conduct 
meaningful temporal autocorrelation analysis. The pC0 2 time series at Hydrostation 
“S” (32°10’N, 64°30’W) and Bermuda Atlantic Time Series Study (BATS) site (31°50’N, 
64°10’W) (Bates et al. 1999) in the Sargasso Sea, did not contain an uninterrupted time 
series longer than a few days. Much of the temporal analysis previously done for OWS 
“P” used, at best, weekly averaged data. For example, Wong and Chan (1991) analyzed 
the seasonal variation of pC0 2 at OWS “P” using monthly and weekly averaged values 
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of pCO s and reported slightly higher values of pC0 2 during summer compared to winter. 
Here the emphasis is on detailed temporal autocorrelation structure. 

The temporal autocorrelation structure for the 11 longest uninterrupted time series 
from OWS “P” are shown in Fig. 5. The 6-hourly reported data has been averaged 
to daily values to remove diurnal effects. All these time series are from the 1970’s with 
the longest segment being 42 days during 1978. Others are as short as 10 days. The 
shorter segments obviously have less confidence associated with the estimated structure. 
Whereas there are instances where the autocorrelation coefficient remains relatively high 
after a 5-day lag, the longer decorrelation scales can really not be associated with any 
particular time of year. The last two structures for the summer of 1978 illustrate this 
point. Whereas for Julian days 171-200 during 1978 the autocorrelation coefficient at 
5-day lag is about 0.6, immediately following for Julian days 202-212, the coefficient is 
near 0 between 2 and 3 day lags. From the data available, it is difficult to assign a general 
structure to the temporal autocorrelation structure of daily pC0 2 , even for a particular 
season. It is possible that longer period averaging, say weekly or monthly, might produce 
more consistent autocorrelation structure, however, the available data can not support 
those calculations. 

Because uninterrupted time series of oceanic pC0 2 are scarce, an attempt is made 
to use other variables as a proxy to estimate the pC0 2 correlation structure, specifically 
temperature because of the great sensitivity of pC0 2 to temperature. Boutin et al. (1999) 
and Cosca et al. (2003) found that some of the pC0 2 variability in the eastern equatorial 
Pacific could be explained by the variability in temperature. Here, the four time series 
with pC0 2 data for at least 30 uninterrupted days is used to investigate the co- variability 
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of pC0 2 with temperature and salinity. The four pC0 2 time series are shown in Fig. 6. 
As in Wong and Chan (1991), most of the values of pC0 2 are between 290 and 340 /xatm. 
The exception occurred during the summer of 1976 when a likely mixing event dropped 
surface pC0 2 values to near 250 /iatm. The temperature and salinity data for this period 
show no unusual behavior. Since no nutrient concentration data are available for this 
period, the low pC0 2 values cannot be attributed unequivocally to biological utilization 
of C0 2 . 

Rather than presenting a simple correlation, we choose to present the coherence spectra 
between the temperature and pC0 2 data. The data is detrended by removing a linear 
best fit line, likely removing most of any local seasonal signal. The estimate of coherence 
between variables is identical to unity except when smoothing is applied to the estimates of 
the auto and cross spectra. For this analysis, a simple Hanning smoother, where a spectral 
component is replaced with a weighted (1-2-1) average of its surrounding components, is 
applied. The coherence spectra for the four time periods between temperature and pC0 2 
are shown in Fig. 7. The value of coherence being 1 at either side of the frequency 
extremes is due to no smoothing of the spectra at end points. The jaggedness in each 
of the coherence spectra is obvious. Whereas the 1976 and 1977 periods might indicate 
greater coherence at lower frequencies, the 1978 sequences do not. Except for the 1977 
series, there is a local peak at 0.2 cycles per day that might indicate that pC0 2 variability 
is related to mixing forced by atmospheric synoptic variability. However, it is difficult to 
say that a clear picture of covariability between pC0 2 and temperature on any particular 
time scale. Although not presented, the results are no more clear for pC0 2 and salinity 
covariability than they are for pC0 2 and temperature covariability at OWS “P”. 
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The seasonal variability of pCC >2 observed in the OWS “P” area is generally less than 
30 /iatm (with an exception of one event in the summer of 1976) as shown in Fig. 6. 
This is about 20% of 140 /iatm that is expected from the 9 °C seasonal SST changes if 
no change in TCO 2 occurred ( d In pCC> 2 /dT = 0.0423 °C -1 ). Takahashi et al. (2002) 
have shown that the temperature effect was largely cancelled by the biological drawdown 
of TCO 2 during spring through summer and by the upwelling of CCVrich deep waters 
during fall through winter. Hence, the C0 2 system in this area depends only weakly on 
the seawater temperature. This accounts for the highly variable coherence between pCOg 
and SST, and the 0.2 cpd peak which may represents the effects of wind-induced deep 
mixing events. 

5. SUMMARY AND CONCLUSIONS 

A global database containing mainly underway measurements of surface water pC0 2 
is used to estimate the spatial autocorrelation structure on a 10° latitude-longitude grid. 
Since the pCOg is a sensitive function of temperature and the total C0 2 concentration 
in seawater, which depends on biological utilization and upwelling of C0 2 -rich subsurface 
waters, its variability reflects both physical as well as biological processes. The PCO 2 
is more densely sampled in space and time in northern hemisphere than the southern 
hemisphere. The longest spatial correlations occur in the equatorial Pacific where the 
autocorrelations are above 0.8 out to distances between 80 — 100km. The other equatorial 
areas also have relatively long correlation scales compared to other areas of the globe. 
The north Atlantic subtropical gyre has spatial correlation scales that tend to be about 
double the correlation scales in the subtropical gyre associated with the Kuroshio. This 
geographic distribution is similar to the distribution of the values of the first baroclinic 
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Rossby radii of deformation as calculated by Emery et al. (1983), indicating that the 
physical properties of the flow field are a likely contributor to factors controlling the 
correlation scales. In the Indian Ocean, the correlation scales in the southern hemisphere 
tend to be longer than those in the northern hemisphere. A possible explanation for 
this variability can again be related to the controlling mechanisms of the flow field. An 
important mode for the southern Indian ocean variability stem from long annual Rossby 
waves that are generated on the eastern side of the basin and propagate westward. The 
northern half of the basin is much more reactive to the strong forcing of the monsoonal 

r' 

flows. The correlation scales tend to be shortest at higher latitudes, again consistent with 
the scales of the deformation radius, but there is an indication that after an initial steep 
drop off of the autocorrelation, there is a slowing of the drop off between 0.8 and 0.6 in 
the northern Atlantic. 

An attempt to estimate the temporal correlation structure is made at the only ge- 
ographic location with relatively long (longer than about 10 days) uninterrupted time 
series of pCOs data, OWS “P”, in the eastern North Pacific. Eleven time series were 
used to calculate the correlation structure with highly variable results. Most of these time 
series indicated a significant drop in autocorrelation (0.5 and lower) after 1-2 days. Three 
of the time series contained a more gentle drop off in the temporal correlation, but after 
5-days, the correlation was only about 0.5. There was no obvious seasonal dependence in 
the temporal autocorrelation structure. Indeed, two of the time series that were separated 
by a single day during the summer of 1978 produced a drastically different autocorrelation 
structure from each other. 
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The coherence of the geographic distribution of the spatial autocorrelation structures 
with physical properties of the flow, such as the internal Rossby radius of deformation, is 
encouraging. As assimilation models are contemplated for pCOz , it is necessary to have 
confidence in an error covariance model for the data being ingested. An autocorrelation 
function provides a good start, and the fact some of the global variability can be related to 
other physical variables gives some hope that any analysis product will be of more use than 
the ingested data alone. However, the inability to estimate a temporal autocorrelation 
structure is troubling. Even at OWS “P” where time series data exists, a stable estimate 
of that correlation structure proved difficult. The situation may not be hopeless. The 
estimates of the spatial correlation structure ignored any time change along a particular 
track, yet those correlations remained high over long distances, which also means over 
a time period of more than a couple days. Although modeling studies may be useful as 
a first guess for estimating pCO s correlation structure, confident estimates will require 
more and targeted measurements over the entire globe. 
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Figure 1. The spatial and temporal distributions of ocean surface pCOs measurements taken 
between 1970 and 2003. The location is marked on the map by a colored marker. The color of 
the marker indicates the pentad during which the observation was taken. 
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white contain no pCOs measurements. 
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Figure 3. (a). Global distribution of lag distance where the spatial autocorrelation exceeds 


0.8, and (b) the lag distance where the spatial autocorrelation exceeds 0.6 
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Figure 4. The estimated average spatial autocorrelation functions for seven specific geographic 


locations (solid lines with circles) and their best fits (dashed lines) assuming an exponential decay 


scale. 
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Figure 6. The four longest time series of daily averaged pCO s from Ocean Weather Station 
(OWS) “P” in the eastern subarctic Pacific (50°N, 145°W). These pC0 2 series are chosen for 
display because the length of each time series is at least 30 days. The year of observation is 
shown above each plot, and the Julian Day is marked on the abcissa. 
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Figure 7. The coherence spectra between temperature and pC0 2 for the four time periods 
shown in Figure 6. The coherence was estimated using a Hanning (1-2-1) filter on the auto- and 
cross-spectral estimates. 
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POPULAR SUMMARY 


Some variables related to climate are not currently observable from space and require 
in situ measurement. The dissolved C02 in the surface water is such a variable and is 
important because it is a prime regulator of how much C02 in the atmosphere is 
exchanged with the ocean. In this manuscript, the natural spatial scales of dissolved C02 
in the global surface ocean are identified. Knowledge of these scales is important as it 
helps determine the requirements for future measurements of C02 variability in the 
ocean. This study clearly shows that the spatial scales of dissolved C02 are intimately 
tied to the scales of the surface ocean circulation. The largest scales occur in the tropical 
Pacific, and shortest scales in polar regions. The estimates of these scales is also 
important for computer models and simulation of C02 variability as these calculations 
provide the first global estimates of the finest spatial scales that models must resolve. 


SIGNIFICANT FINDING 

Although there have been a handful of studies documenting the observed scales of 
ocean surface C02, this is the first study to present results of ocean C02 scales for the 
entire globe. 



